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I. INTRODUCTION 



The E821 experiment at the Brookhaven National Laboratory has recently measured the anomalous magnetic 
moment of the muon, a M = (<?„ — 2)/2, with the final result 

a e * p = 11 659 208.0(6.3) ■ 10~ 10 . (1) 

This unprecedented accuracy, with yet better precision expected in the planned experiments at BNL, JPARC, and 
FNAL [2|, maintains the live interest in obtaining more accurate theoretical predictions for a M within the standard 
model, for reviews see, e.g., 0, 0) IS Q- The challenge is to obtain the theoretical uncertainty lower than the 
uncertainties for the nearest-future experiments, which will supply a powerful test for possible effects of contributions 
from the New Physics. 

In the standard model the electromagnetic (QED), electroweak (EW), and hadronic (Had) effects contribute to a M . 
Taking into account the recent progress with the QED calculations and the latest result for a e one obtains [7| 

a QED = n 658 471.809(0.016) • lO" 10 . (2) 

The electroweak contribution to a M is also known accurately [f| , 

a™= 15.4(0.2) • 10~ 10 . (3) 

The main source of the theoretical uncertainty is the hadronic contribution. There are three types of the leading 
hadronic contributions: the vacuum polarization, its next-to-leading folding with the QED and EW sectors, and the 
light-by- light (LbL) scattering process (Fig. |lj) . A recent phenomenological reanalysis of the contribution of the full 
hadron vacuum polarization insertion into the electromagnetic vertex of the muon 0] gives 

a Had,LO = 69(I9 (44) . 10 -io ( 4 ) 
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The most recent estimate of the higher-order (HO) hadronic contributions performed in flQj provides the result 



a 



Had,HO_ 
I' 

For the LbL contributions there are several model-dependent estimates: 



= -9.8(0.1) • 1CT 10 . (5) 



Had.LbL = 8.3(3.2) • 10" 10 [11,12], (6) 
Had,LbL = 8 .9(1.7) . 10" 10 [13], (7) 
Had,LbL _ i3. 6 (2.5) • 10" 10 [14]. (8) 
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(a) (b) 

FIG. 1: (Color online) (a) Hadronic light-by-light contribution to a M . The bottom line is the muon, the wavy lines are the 
photons, and the circle depicts the hadronic part, (b) The meson pole contribution to a M . The circles represent the virtual 
meson to 7*7* transition form factors. 



It is clear that the leading-order (LO) hadronic contribution (HJ) dominates in the absolute value. However, the 
theoretical error introduced by the LbL process |6])-([8]) is of the same order as that for the leading term. Moreover, 
the error in ([4]) is phenomenologically well controlled and may be improved by a factor of 2 or so if more experimental 
data on e + e _ (t) — > hadrons should appear [9J. The precision for a^ ad,HO is quite enough for the nearest-future 
experiments. Unfortunately, the hadronic LbL contribution cannot be related to any other observable and hence we 
must rel y o n purely theoretical model framework in order to estimate it. Predictions based on chiral models ([6]), d7]) 
(also in [jjl) are compatible with one another and are much lower than the results obtained in (JSJ) (also in (lq, fl7j|). 
Furthermore, the errors given in ©-([5]) are difficult to estimate and it is not easy to improve their quality. These 
errors should include the model dependence, the dependence on model parameters, and the quality of the model 
assumptions. The last point, difficult to assess, is most important for our further discussion presented in this pa per . 

We digress that at present the calculation of the LbL contribution within lattice QCD is a quite difficult task [l8| , 
hence we cannot use lattices yet as a reliable source of information for the considered problem. 

The present work is devoted to the calculation of the pion-pole contribution of the hadronic light-by-light scattering 
to within the nonlocal chiral quark model (NxQM) [lij [2(j motivated by the instanton model of the QCD 
vacuum [U [22, HH. In [2(| the vector and axial- vector correlators were calculated and then these results were 
applied to compute the hadronic vacuu m p olarization contribution a^ ad ' LO [24J. Later on the three-point PVV and 
VVA correlators were analyzed in [25|, [III [13] and the contribution of the 77* Z* vertex to was estimated in 
The present calculation of the LbL contribution exhibits a few important improvements compared to the previous 
calculations in effective quark models. First of all, in the pion-pole contribution, which dominates the hadronic LbL 
part, we take into account the full kinematic dependence of the pion-photon form factors, including their dependence 
on the pion virtuality. The inclusion of this dependence diminishes the results for the contribution to a M by about 
15% compared to the calculation with no dependence on the pion virtuality. Importantly, our approach is consistent 
with the low-energy theorems and with the QCD constraints. In particular, we demonstrate how the QCD constraint 
considered by Melnikov and Vainshtein |l4[ is satisfied at high photon momenta within N^QM, with the crucial role of 
the box diagram. We also show that in general the pion-pole contribution exhibits no enhancement as it was assumed 
in the model [3| for the pion exchange in the LbL amplitude. 

The next Section contains the definitions of the LbL amplitude, its pion-pole contribution, and the QCD constraint 
for the amplitude discussed in [3]. Sections III and IV describe the nonlocal chiral quark model based on instanton- 
like dynamics (NxQM), introduces conserved vector and axial- vector currents, as well as the relevant PVV amplitude. 
This part, developed in previous works, is included for the completeness of the paper. Numerical results for the pion 
pole contribution to LbL and the comparison to other calculations are shown in Sections V and VI. In Section VII 
the Melnikov- Vainshtein constraint is proven within NxQM. Conclusions are given in the last Section. 
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II. LIGHT-BY-LIGHT AMPLITUDE 

The amplitude for the light- by- light scattering is defined as 1 

M = a 2 N c Tr [Q 4 ] A = a 2 N c Tr [Q 4 ] A^^ge^ ef ef f* 

= -e 3 f dWt/e-^-^ef tftf < 0|T (x) j M2 (y) j M3 (0)} | 7 >, (9) 



where <^ and ej are the momenta and the polarization vectors of photons, is the hadronic electromagnetic current 
defined explicitly below within the nonlocal chiral quark model, and Q is the quark charge operator. The photon 
momenta are taken to be incoming, with ^2 Qi — 0. One of the photons represents the external magnetic field and can 
be regarded as a real photon with a vanishingly small momentum, q^. Due to the gauge invariance the light- by- light 
scattering amplitude is proportional to the field strength tensor of the soft photon, f lS = qle.\ — qfej. Neglecting the 
quadratic and higher powers of (74, the tensor amplitude A^^^^s may be considered as a function of the photon 
virtualities qf, z = 1,2,3 under the condition qi + q2 + 93 = 0. 

In general the LbL amplitude is a rather complicated object to analyze. However, it is possible to generate different 
hierarchies for the components of the amplitude. The dominant contribution comes from the pole in the pseudoscalar 
channel, which is enhanced by the very small value of the pion mass (Fig. [TJ)). It is also leading in terms of the 
large--ZV c counting, where N c is the number of colors. Moreover, the leading-7V c contributions from other mesonic 
channels are suppressed by much larger meson masses. 

It can be shown, see [lj, [HI, [28[ , that the leading contribution from the neutral pseudoscalar meson exchange to 
a M is given by 



a LbL;x" = _ e 6 



d i q l f d 4 q 2 



' / (2nyj (2tt)4 qlqlql[{j> + qxy-m?][{p-q 2 y-m?] 

F-** TT (ihQi^D-n-r, — °n T , 2\\ T KTT (&<&,0) Ti(qi,q 2 ;p) 

9-xq - <^P J pp \12)) 

+ ■ ? >37'7*(93;9l)92)-2-Ti — °n r 7 2^ ^1*1* (ih (li + I2) 2 , 0)T 2 (q 1: q 2 ;p) 



(10) 



where m denotes the muon mass (p 2 — m 2 ) and the kinematic factors are [HI 

Ti(qi,q 2 ;p) = y (p ■ qi) (p ■ q 2 ) (<?i • ftO - y {p ■ q 2 ) 2 qf 
8 16 

~ 3 (p • 9i) {qi ■ 92) ql + Hp ■ 92) q\ q\ - y (p ■ 92) (<?i • <?2) 2 

+ y m 2 q\ q\ -ym 2 (qi • q 2 f , (11) 
T 2 (qi,q 2 ;p) = y (p ■ qi) (p • q 2 ) (qi ■ q 2 ) - y {p ■ qi) 2 ql 
+ \(p- <?i) (<?i • 92) ql + I • qi) q\ q\ 

+ \rr?q\q\ - \m 2 { qi ■ q 2 ) 2 . (12) 

The PVV amplitude ^■•7*7* (<?f; q 2 , q 2 ) for the virtual pion and photons, as well as the meson propagator in the 

pseudoscalar channel, (l — GpJ pp (q 2 )) , will be discussed in Sect. IIVI 

Recently, an important constraint was introduced by Melnikov and Vainshtein (MV) [HI , who argued with the help 
of the operator product expansion (OPE) that in the specific kinematic limit, q\ ps q\ = q 2 ^> q 2 , the amplitude has 
the asymptotic form (see Fig. [2]) 

A^ 3l sP S = ^ llJi2Sp q s W(a) {™ ( L%i) tiQzUz + ™t%1) k~<lkU^) } + ■••, (13) 

5 a=3,8,0 



1 In the following definitions we follow closely the notation used in Il4l . 
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where f a p = ^e al 3 7 sf lS , Q = (li — 92) /2 and no hierarchy between q 2 and Aq CD is assumed. The weights are 
defined as 

2 

l2l 



Tr [\ a Q 

w (a) = \ J_ / 14 N 

Tr[A2]Tr[g4] 

W (V = I W (») = J_ ^(O) = H 
4' 12' 3' 

The invariant functions w^\q 2 ) and w^\q 2 ) are related to the triangle amplitude that involves the axial current 
ji^p and two electromagnetic currents, one with momentum q and the other one (the external magnetic field) with 
the vanishing momentum (see the relevant part of Fig. (2b) , 

T$= * (0| / d 4 ze^T{j^(z) j M (0)}| 7 >. (15) 

It is shown in f2^ | that t!$ can be written as 

+4V)(-'? 2 /mp+9^ ,T /^-9^ ,T /^) } • (16) 

The first (second) amplitude is related to the longitudinal (transverse) part of the axial current, respectively. In the 
chiral limit wl is fixed by the axial anomaly and the non-renormalization theorem predicts [29j 

-,„( 3 > 8 ) f„2\ 



2/q*. (17) 

QCD 



In the hard limit q 2 3> Aq CD the perturbative QCD predicts the asymptotic form (a = 3, 8, 0) 



4° (l 2 ) = (l 2 ) = 1/9 2 (18) 

and (|13| transforms into 

AiM2M37<5/ 7<5 = -^ e MiM2<5p9 A { 2 9393/^3 +(-93/^+93^93/^-9393/^3)} + ■■■■ (19) 
93 9 x ' ' 

It was demonstrated in [2(| [27J that the results (|16til8|) are satisfied within the non-perturbative N%QM approach. 

In the next sections we discuss the basic elements of the NxQM. Then within this model we calculate the leading 
pion-pole contribution to a^, demonstrate how the Melnikov-Vainshtein asymptotics is realized, and finally discuss 
on the role of the triangle functions u>l(t) for a^. 
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III. NON-LOCAL CHIRAL QUARK MODEL 

To study the LbL contribution to a M one can use the framework of the effective field model of QCD. In the low- 
momenta domain the effects of the non-perturbative structure of the QCD vacuum become dominant. Since the 
invention of the QCD sum rule method based on the use of the standard OPE, it is common to parameterize the 
non-perturbative properties of the QCD vacuum in terms of infinite towers of the vacuum expectation values of the 
quark and gluon operators. From this point of view the nonlocal properties of the QCD vacuum result from the 
partial resummation of the infinite series of power corrections, related to vacuum averages of the quark an d g luon 
operators of growing dimension, and may be parametrized in terms of the nonlocal vacuum condensates [30l. I3lll32l]. 
This construction leads effectively to nonlocal modifications of the propagators and effective vertices of the quark and 
gluon fields at small momenta. 

An adequate model embodied in this general picture is the instanton model of the QCD vacuu m |2lL I22L I23I ] which 
describes non-perturbative nonlocal interactions in terms of the effective action (for a review see 33]). Spontaneous 
breaking of the chiral symmetry and the dynamical generation of a momentum-dependent quark mass are naturally 
explained within the instanton model. The non-singlet and s ing let vector and axial-vector current-current correlators 
and the vector Adler function have been calculated in [2(J Hj| m the framework of the effective chiral model with 
instanton-like nonlocal quark-quark interaction [li^ . |20| . It was shown, in particular, that the Weinberg sum rules 
which are sensitive to the low- as well to the high-momentum dynamics are satisfied [20I l3~il . In the same model 
the pion transition form factor normalized by the axial anomaly has been considered in [25| for arbitrary photon 
virtualities. 

We start with the nonlocal chirally invariant action 2 which describes the interaction of soft quark fields [36j . 



S = J d 4 x qj(x) \%>fD M -m f ] qi {x) + (20) 

+ ^G P J d 4 X j Y[d 4 x n f(x n )\Q(X- Xl ,X) ■ 

■ T P Q(X, X + x 3 )Q(X - X2,X)T P Q(X, X + x 4 )] , 

where m/ is the current quark mass, — — iV^ (x) — ij^A^ (x) is the covariant derivative, and the matrix product 
rpigTp = (1 (g> 1 + i-fsT a (g) i-fsT a ) provides the spin-flavor structure of the interaction. In Eq. (f2"0")) qj = (m, d) denotes 
the flavor doublet field of dynamically generated quarks, Gp is the four-quark coupling constant, and t" are the Pauli 
isospin matrices. The separable nonlocal kernel of the interaction described in terms of the form factors f{x) is 
motivated by the instanton model of the QCD vacuum, where the function f(x) may be evaluated. In order to make 
the nonlocal action gauge- invariant with respect to the external vector and axial gauge fields V^{x) and A®(x), we 
employ in (12T))) the delocalized quark field, Q(x), defined with the help of the Schwinger gauge phase factors, 



Q(x,y) =Pexp|* J d Zfl [V£(z) + l5 A«(z)]T a j qi (y), 

Q(x,y)=Q\x,yh°. (21) 

Here P stands for the operator ordering along the integration path, with y denoting the position of the quark and x 
being an arbitrary reference point. 

The dressed quark propagator, S(p), is found to be 3 

S- 1 (p) = ip-M(p 2 ), (22) 

with the momentum-dependent quark mass defined as the solution of the gap equation 

Mtf) = m f + 4G P N f N c f(j>*) J |* f^) fc2 ^| fc2) . (23) 



2 In the present work we do not consider extensions of the model to include explicit vector mesons and the effects of the strange quark 
mass. Such extensions of the nonlocal model were considered in [53 . |35|| and will be used in complete calculations elsewhere. 

3 From here on the Euclidean metric for the momenta is used. 
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The formal solution is expressed as [33] 

M{p 2 ) = mj + (M q - m f )f 2 (p 2 ), (24) 

with the constant M q = M(0) determined dynamically from Eq. (123p . The momentum-dependent function f(p) is the 
normalized four-dimensional Fourier transform of f(x) given in the coordinate representation. Within the instanton 
model / is related to the zero mode solution for the massless fermion in the field of an instanton. In [26( it was shown 
that in the momentum space the nonlocal function / (p) for large momenta must decrease faster than any inverse 
power of p 2 , e.g., like an exponential, and that this choice corresponds to calculations in the axial gauge for the 
quark effective field. In order to take these effects into account and to make numerics simpler we use for the nonlocal 
function the Gaussian form, 

f{p) = exp (~p 2 /A 2 ) , (25) 

where the parameter A characterizes the nonlocality size of the gluon vacuum fluctuations and is proportional to the 
inverse average size of the instanton in the QCD vacuum. 

An important property of the dynamical mass (|23[) is that at low virtualities its value is close to the constituent 
mass, while at large virtualities it goes to the current mass value. This property is crucial in obtaining the correct, 
consistent with the OPE, QCD behavior of different correlators at large momentum transfers [2(| HH HH • The nonlocal 
chiral quark model can be viewed as an approximation to the large- N c QCD, where the only (effective) interaction 
terms, retained after integrating out the high-frequency modes of the quark and gluon fields down to the nonlocality 
scale A where the spontaneous chiral symmetry breaking occurs, are those which can be cast in the form of the 
four- fermion operators (|20|). The parameters of the model are then the nonlocality scale A, the four- fermion coupling 
constant Gp, and the current quark masses m/. 

The quark-antiquark scattering matrix in the pseudoscalar channel is found from the Bethe-Salpeter equation as 

f '<' 2 > - i-cwr (26) 

with the polarization operator in the pseudoscalar channel equal to 

Jpp(q 2 ) = [ t^V 2 (*0 f (k + q) Tr [S(k) 75 S (k + a) 75 ] . (27) 

J (27T) 

The position of the pion state is determined as the pole of the scattering matrix 

det(l - GpJ PP (q 2 ))\ q2= _ m , = 0. (28) 

The quark-pion vertex identified with the residue of the scattering matrix is (k' = k + q) 

r*(k,k')= gTrqq t l5 T a f(k)f(k') (29) 



with the quark-pion coupling equal to 



_ 2 _ dJpp (q 2 ) 



) 71 q 



dq 2 



(30) 



where is the physical mass of the pion. In the chiral limit the quark-pion coupling, g vq , and the pion decay 
constant, /„■, are connected by the Goldberger-Treiman relation, g n — M q / f^, which is verified to be valid in N%QM, 
as requested by the chiral symmetry. 

IV. CONSERVED VECTOR AND AXIAL- VECTOR CURRENTS 

The vector vertex following from the model (|2"0|) is 

T^(k,k')= Q\ lfi + {k + k%M^{k,k')], (31) 
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where Q is the diagonal matrix of the quark electric charges, AfW(fc, k') is the finite-difference derivative of the 
dynamical quark mass, k (k 1 ) is the incoming (outgoing) momentum of the quark, and q = k' — k is the momentum 
corresponding to the current. The finite-difference derivative of an arbitrary function F is defined as 



(32) 



The full axial vertex corresponding to the conserved axial-vector current is obtained after the resummation of quark- 
loop chain that results in the form of the term proportional to the pion propagator [19] 



rf(k,k>) = [ 1 ^ 5 + Arl(k 1 k')]x a , 



AT* (*,*') =2 75 %/(fc)/(fc' 



Jap (0) - 



rrifGpJp (q 2 ) 
1 - GpJpp (q 2 ) 



+ -f 5 (k + k')^J AP (0) 



(f(k')-f(k))' 
k' 2 - k 2 



(33) 



where we have introduced the notation 



Mi 2 ) = J 



d A k 



f(k)f(k + q)Tr [S(kh 5 S(k + q) l5 ] 



(34) 



J AP (°- 2 )=4N c N f J £J ™W y/M(l + q)M(l). 



In the chiral limit of to/ = one gets instead 



Alt(fc, k') = 75 % [M(k') + M (k)} + 75 



it 



(k + k') f 



k' 2 - k 2 



M q (f(k')-f(k)Y 



The axial-vector vertex (f3"3"|) has a pole at 



Q 2 = -ml = to/ (qq) /ft 



(35) 



(36) 



(37) 



where the Goldberger-Treiman relation and the definition of the quark condensate has been used. The pole is related 
to the denominator 1 — GpJpp (q 2 ) in Eq. (|33[) . while q 2 in denominator is compensated by zero from the square 
bracket in the limit q 2 — > 0. This compensation follows from expansion of J(q 2 ) functions near zero momentum 



— G 



Jpp(q 2 ) 
J A p(q 2 = 0) 

It follows from this expansion that near pole position, q 2 

G P 



to/ (qq) M- 2 - q 2 g- 2 + O (q 4 ) , 
M q , J P (q i = 0) = (qq)M-\ 



-to 2 , one has for the factor in (flT))) 



1 



(38) 



gl „(1- GpJ pp (q 2 )) q 2 + 



(39) 



However, at large q 2 this factor goes to the constant Gpg~ 2 instead of decay like q~ 2 as it would be predicted by the 
simple pole approximation. In the chiral limit m/ — the second structure in square brackets in Eq. (|33[) disappears 
and the pole moves to zero. It should be stressed that the form of the interaction vertices (|31I33[) is consistent with 
the chiral Ward-Takahashi identities. 

Within NxQM the full singlet axial-vector vertex including local and nonlocal pieces is given by [2(| 



r°' 5 (M') =7p75 + Ar°' 5 (fc,n 



(40) 



Ar^(M') 



75 



(k + k') p M q 



(f(k')-f(k)) 2 
k' 2 - k 2 



C ^2M q f{k')f{k) G ^ l - GpJpp(q2 



G P 1 - G p Jp P {q 2 ) 



where G P is the four-quark coupling for the singlet channel. As it follows from expansion |38|) the singlet current ([40]) 
does not contain massless pole due to presence of the Ua (1) anomaly. Instead, the singlet current develops a pole at 
the rj'— meson mass 



l-G^Jpp( g 2 = -TO 2 ,) = 0, 



(41) 
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thus solving the I/a(1) problem. 

The parameters of the model are fixed in a way typical for effective low-energy quark models. One usually fits the 
pion decay constant, /„■, and the pion mass (]28[) to their experimental values. In the chiral limit the decay constant 
reduces to fo^ — 86 MeV (38|. In NxQM /o jjr is expressed as [23j 



/o,7r 



4tt 2 



M 2 (fc 2 ) - k 2 M(k 2 )M'(k 2 ) + k 4 M\k 2 ) 2 
(fc 2 + M 2 (k 2 )f 



dk 2 k 



(42) 



where the primes mean the derivatives with respect to k 2 , i.e. M'(k 2 ) — dM (k 2 ) / dk 2 , etc. The described fitting 
procedure introduces two relations among the three model parameters. We use the following values of the model 
parameters fixed in Ref. [zjj : 



M q = 0.24 GeV, A P = 1.11 GeV, m/ = 8 MeV. 
To test the sensitivity of the results on the choice of parameters we also use the set [If 

M q = 0.35 GeV, A P = 1.2 GeV, m f = 12 MeV. 



(43) 



(44) 



Both sets reproduce the low-energy observables with an acceptable accuracy of the order of 10%. The finite current 
quark mass in (|43|) and (|44| is fixed by the physical value of pion mass. It also increases the pion decay constant to 
its experimental value, f v — 92 MeV, within a 2 MeV accuracy. 







+ 




FIG. 3: (Color online) The diagrammatic representation of the LbL scattering amplitude within the effective four-quark model 
in the leading 1/N C approximation. It consists of the box diagram plus the iteration of the four-quark interaction term via 
the quark loop. The iterative terms sum up into a two-point meson correlator in the given channel (pseudoscalar, scalar, or 
axial- vector) . All quark lines and the vertices are dressed, as calculated within NxQM. There are also a number of contact 
terms inherent to nonlocal models, not shown in the figure. 



The triangle PVV amplitude corresponding to the process 7*7* — > 7Tq 

A (7* ei) 7* (Q2, e 2 ) -> tt° (p)) = -ie 2 e^ vpa ^qlqlT^ aTT {(g\ q 2 l7 gg), (45) 

has been constructed and discussed in [25j 4 . In (|45f and qi (i = 1,2) are the photon polarization vectors and 
momenta, while qg, = q\ + q2 is the pion momentum. In NxQM one finds the contribution of the triangle diagram to 
the invariant amplitude as 



A (7*72 



-ie -^—1 r 



Q 2 t 3 



(9i 



q 2 ;ei <-> e 2 ] 



ef e 2 



d 4 fc 



Tr[Yl (fc+, 5(fc_)r M (fc_, k 3 )S (fc 3 ) r„(fc 3j k+)S(k+)} (46) 



4 In |25| the amplitude T v 



is denoted as 
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where q — qx — q 2 , k± = k ± q 3 /2, k 3 — k — q/2. In the chiral limit (g| = m 2 . = 0) with both photons real [q 2 = 0) 
one finds the normalization by the axial anomaly (within the nonlocal models it was proven in [35j ) 

^, 7 , r (0;0,0) = I ^r-. (47) 

The functions defined in (|26[) . ([27)l . (|4"5)l . are used in the evaluation of p0[) . Within the effective model the 
light-by-light amlitudc is given by diagrams of Fig. [31 



V. NUMERICAL RESULTS FOR THE PION POLE CONTRIBUTION TO LBL 



We first investigate the role of the off-shellness in the pion transition form factor JT 7r * 7 , 7 » (q|; q\^ g|), which is the 
important element considered in this work. In Fig. 0] we show how T^*^^* (gf; q 2 , q 2 ) I Fit'-i*^' (0; 0, 0) depends on q 3 
for fixed values of q\ and q 2 . For all cases of q\ and q 2 we note a strong dependence on q 3 , which for large values of 
this momentum becomes exponential, ~ exp(— const • eg). This strong dependence on the off-shellness, related to the 
Gaussian form of the nonlocal quark-pion vertex (f!?9"]) , could in principle lead to a significant reduction of the pion 
pole contribution to LbL. As we shall see below, this is not the case, since the effective support for the integrand 
of a^ bL,7r ° is localized at moderate values of q 2 where the form factor suppression is not so strong. Furthermore, 
the factor (l — GpJ pp (q 2 )) 1 in Eq. (TTUJ) at large q 2 goes to one, but not to q 3 2 as in the case of the single-pole 
approximation (q 2 — mj) . 




FIG. 4: (Color online) The dependence of the pion transition form factor T^*-,*-,* (<j§ ; q\ , gf ) / 'J r 7r * 7 » 7 » (0; 0, 0) on the off-shellness 
of the pion, |g 3 |, for q\ = q\ = (top curve), q\ = (700 MeV) 2 , q\ = (middle curve), and for q\ = q\ = (700 MeV) 2 /2 
(bottom curve). We note a strong dependence on \qz\ (all momenta are Euclidean). 



To obtain a^ bL:7r °, one has to convolute the three-point vertex functions discussed above in the two-loop integral of 
Eq. (jlOp . Since the multidimensional integrand is a regular function, this task can be accomplished straightforwardly 
with the Monte Carlo integration technique. We use the routine VEGAS H^|, which implements the adaptive Monte 
Carlo algorithm. The eight-dimensional 5 integral in (|I0p . carried over the the two muon loop momenta, has been 
done numerically in the Euclidean space. 

A convenient and important trick, allowing to avoid the numerical cut-offs in the momentum integration, is to map 
the momentum variables via a conformal transformation into the range [0, 1). Explicitly, we use 

C -TZ7> ( 48 ) 



P 

with pi denoting the Euclidean momentum and a = 1 GeV 2 . We have tested the accuracy of our code by reproducing 
very accurately various model results listed in Ref. [l5[. 

The result of the calculation in N%QM with parameters (|43|) is 



5 Through the use of the rotational symmetry one may reduce the integral down to five dimensions, but from the view point of the Monte 
Carlo integration this is irrelevant for the performance. 
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while with parameters (|44|) we get 

a LbL^o = 6 . 68 . i - 10 . (50) 

Similarly to other models, our result is dominated by the term proportional to the T\ structure in Eq. (I10p . which 
yields 97% of the total of (129150")) . 

In order to illustrate the convergence of the result with the increasing upper limit of the momentum integration, 
here denoted by A, we show in Fig. [5] the relative yield of the full result for ap bL,7r ° with model parameters (|44|). 
plotted as a function of the Euclidean momentum cut-off, A (\pt \ < A). For A — > oo we obtain, by definition, the full 
result. For A = 1 GeV about 90% of the full result is obtained. Thus the result is dominated by the soft physics, as 
requested of the effective low-energy model. These findings are in agreement with that obtained in [Til, |4(| . 



Relative yield 

1.00 




-■- A [GeV] 

3.0 



FIG. 5: (Color online) The relative yield of the full result for a M ,1r ° for model parameters (|44|) . plotted as a function of the 
Euclidean momentum cut-off, A. 



Finally, we note that in the extended version of the considered model the vector and axial- vector couplings appear. 
However, as it was discussed in [25[ within the nonlocal model, the inclusion of these couplings leads to very minor 
changes of the pion transition form factor. In particular, it does not change the normalization of the form factor at 
zero momentum and contributes small corrections to the leading-order asymptotics. The latter is due to the relatively 
large mass of the rho meson. 



VI. COMPARISON TO OTHER APPROACHES 



In this section we compare our calculation of the pion-pole light-by-light contribution to a M with the results obtained 
in earlier calculations [ill, EE EH 



_0 

a 



.jr°,LbL = 5. 6 .10" 10 [11], (51) 

a £°,LbL =5 . 6 .io- 10 [13], (52) 

a£°< LbL = 5.8(1.0) -10- 10 [15]. (53) 

These calculations are based on the usage of different parameterizations of the 7r7*7* vertex satisfying the CLEO data 
on the pion transition form factor 7r — > 77* and OPE constraints. Thus, in [151 ] the phenomenological form factor 
from the generalized vector-meson dominance (VMD) was used (the Knecht-Nyffeler (KN) model), 



pgVMD, .x _/* (« + t) st - h 2 st + h 5 (s + t) + M v M v h 7 
*Tf [S ' ' 3 (M v + 8 ) {My + t) (M Vi + s) (M Vi + 1) ' 

with the parameters My = 769 MeV, M Vl = 1465 MeV, h 5 = 6.93 GeV 4 , h 7 = N c / (^ 2 fl). The error in © 
is due to the uncertainty of the parameter h 2 taken from the interval h 2 £ [10, —10] GeV 2 . However, the function 
(f54j) depends on two kinematic variables instead of three required by (flU)) . We test the effect of this approximation. 
Namely, in (|46|) we have removed the dependence on the pion virtuality and have taken the meson propagator in (I10|) 
in the single-pole approximation 1/ [p 2 — m 2 ^. In the NxQM we obtain numerically for the parameter set (|43[) 

a LbL,^ = 7 33 . 10 -io ^ no off _ shellness effects). (55) 
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This result should be compared to the full calculation yielding (|49f . Thus, as already stated, the inclusion of the 
dependence on the pion virtuality reduces the result by about 15%. A similar statement that the influence of the pion 
off-shellness has only a small effect was presented in [ijjj. In that work different form factors, including the point-like 
pion coupling and the form factors used in the extended Nambu-Jona-Lasinio (ENJL) model of [HI 03) were tested. 
The stability of the results was demonstrated. 

Recently, there was an attempt by Jegerlehner and Nyffeler [4l| to improve the KN model by taking into account 
the full kinematic dependence of the pion-photon vertex through the use of the parametrization for i 7 '^» 7 » 7 . (s,i, u) 

suggested in (42J. This yields for the sum of the tt, 77, 77' contributions a)j bh ' 7To ' T, ' r ' — 9.66 ± 4.5 • 1CP 10 . As a result 
an enhancement of the contribution in comparison with KN result o^.'o.u.u _ 8 3 ± l 2 . i -w 

was obtained. We 

believe that this is a result of neglecting the dressing effects in the parametrization constructed in [42|. In particular, 
this parametrization reproduces the short-distance behavior of the three-point correlator of the pseudoscalar and 
vector currents [42|. For these currents the local 75 and 7 M couplings were used, correspondingly. Because of the local 
character of these vertices the correlator has a slow power-like behaviour at large momenta. However, this correlator 
does not correspond to the triangle quark diagram connecting the virtual pion with the photons. The pion has its 
own hadronic form factor, which means that the 75 vertex is dressed, as it corresponds to the nonlocal form factor 
(|2"9"| (see the similar results in the Schwinger-Dyson approach [43j]), As a result the large- momentum properties of 
the calculation change drastically. 

Similar problems are encountered in calculations based on the usage of the constant constituent quark mass (the 
local models) (l6l . [l7| . Typically, these models predict larger numbers for the hadronic contribution due to incorrect 
form factors, artificially enhanced with respect to the OPE asymptotics. 

We have to note that our approach is closest to the formalism of the work pjj] (ENJL model) in the sense of the 
use of techniques of the effective models. However, the physical grounds are quite different. The ENJL model is in 
fact a generalization of the vector meson dominance model, while N^QM is based on the fundamental property of 
the QCD vacuum, namely the nonlocality of the vacuum fluctuations of the gluon field . Unfortunately, the region of 
applicability of the ENJL model is restricted to momenta much lower than 1 GeV. For example, we are not able to 
prove the properties of the triangle diagram (|16H18[) by using the VVA functions found in the ENJL model [12] ■ The 
inconsistency of the Adler function found within ENJL with the OPE was noted earlier in [2J] . The high-momentum 
region was then modeled in (l2| by different parameterizations, which may be quite risky. One of us (AD) has already 
pointed out in (24| that a possible reason for the failure of ENJL at large momenta is due to the fact that this model 
is based on the single resonance approximation, while nonlocalities, inherent in N^QM, effectively take into account 
an infinite number of resonances according to the quark-hadron duality principle. One of the results that is observed 
in [12] is that there are large cancellations of different contributions. This might be a model-dependent statement 
that has to be checked in future complete calculations within the nonlocal chiral quark model with instanton-likc 
interactions. 



VII. THE MELNIKOV-VAINSHTEIN CONSTRAINT IN N X QM 

The sole contribution of the pion pole discussed in the previous sections and also used in [ll[ and [l5j] does not 
satisfy the Melnikov-Vainshtein (MV) asymptotics (fT9|) discussed in Sect. ITT] In the present section we show within the 
NxQM model that it is the quark box diagram of the light-by-light scattering that provides the correct asymptotics 
(Fig. To this end we first recall the properties of the non-diagonal correlator of the vector current and the 

axial- vector current in the external electromagnetic field (VAU) (fT5|) . as derived in [2(| H3| within N%QM. In this 
model the VAV correlator is defined by (Fig. [5]) 



T; 3p = -2eN c Tr 



d A k 
(2^ 



Tr [r M3 {k + q 3 ,k)S(k + q 3 ) 



■T a / (k + q 3 ,k-q 4 )S(k~ q 4 ) T v (k, k - q 4 ) S (k) 



(56) 



where for definiteness we consider 54 as a soft momentum. The quark propagator, the vector vertex, and the nonsinglet 
and singlet axial- vector vertices are given in (|22|) . (|3Tj) . (f33|) . and (f40|) . respectively. The structure of the vector vertices 
guarantees that the amplitude is transverse with respect to the vector indices 



T a o M3 — D 
1 ii 3 p% — u ' 



(57) 



12 




FIG. 6: (Color online) Diagrammatic repre- 
sentation of the triangle diagram in the in- 
stanton model with dressed quark lines and 
full quark-current vertices 



1,0 

Q\(0> 

0,8 



Q(GeV) 



FIG. 7: (Color online) The normalized nonsin- 
glet invariant function wl constrained by the 
Adler-Bell-Jackiw anomaly following from the 
triangle diagram. The dashed line denotes the 
local part and the dash-dotted line the nonlo- 
cal part. 



while the Lorentz structure of the amplitude is given by ([16]) . It is convenient to separate the correlator into three 
pieces: 



P3P 



a^a,Polc 

P3P 



a^a,Rc st 
P3P 



-2eN r Tr 



-2eN c Tr 



-i2 \a 



d 4 fc 

(2^ 
d 4 k 



Tr [Jh3 S ( k + 93) IplsS (k - q 4 ) j u S (k)} 



(58) 



r M3 (* + Q3, k)S{k + q 3 ) AT 5 p (k + 93, k - q 4 ) S(k- q 4 ) T u (k, k - q 4 ) S 



rpa 
M3P 



rpa,hoc rpa,Po\e 



(60) 



In Fig. [7J it is shown (for details sec 26]) how the different terms saturate the anomalous nonsinglet amplitude wl 
(fT7|) . The part T^^° c , where all vertices are local, saturates the anomaly at large momenta, which is just the property 
of the quark triangle diagram (dashed line in Fig. [7]) 



rpa,Loc I 

M3P 14 soft 

93 >1 GeV 



rpa 



(61) 



This result (|6"Tj) is independent of the channel and the relations (fT8| become true (see [2^, [23] for details). 

It is the pole term (I59[) which dominates the low-momentum behavior of the amplitude sensitive to the flavor 
structure of the axial current. In the nonsinglet channel at zero momentum the anomaly is saturated by the massless 
pion pole contribution (I59| (dash-dotted line in Fig. [Jj). The remaining part T^ cst , which vanishes at zero and 
infinite momenta and is numerically small everywhere, accomplishes the exact saturation to the correct value at 
all momenta. On the other hand, in the singlet case due to different structure of the singlet current (|40[) the pole 
contribution (|59[) is almost completely compensated and the local part (|58[) dominates the amplitude at all momenta 

53- 

With this background in mind let us now consider the LbL scattering amplitude coming from the box diagram with 
dynamical quarks (q 4 is a soft momentum) 



yiBox jyy6 P4 

PlM2P37'5'' — i 



J ^Tr [5(A;)r M4 (k,k + q 4 )S 



(k + q 4 ) T M1 (k + q 4 ,k + q 4 + g x ) S (k + q 4 + q x ) (62) 



r M2 [k + q 4 + qi, k - q 3 ) S (k - q 3 ) r„ 3 (k - q 3 , k)] 



The kinematic limit considered by Melnikov and Vainshtein (q 2 ~ Q2 = 9 2 ^ I3) 1S very 

pion transition form factor with highly virtual photons analyzed within NxQM in (25L l44| . In this limit one has 



q\ = q 2 3> q 3 ) is very similar to the case of the 



^Box ryS 



2tt i f d k 

-^-£pip 2 5p^ e 4 4 J ~^i Tr \- S r P4 ( fc , k + q 4 )S(k + q 4 ) 7^756" (k - q 3 ) (k - q 3 ,k)] - 



., (63) 
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where the higher power corrections are denoted by dots. The result (|63l) is proportional to the correlator T^ oc which, 
as discussed above, saturates the full triangle amplitude at q 2 3> 1 GeV 2 . Thus, within NxQM it is the box diagram 



which saturates the MV large-g 2 asymptotics (fP9|) at q 2 ^> 1 GeV 2 



A 2 
QCD • 



At the same time at q§ <C 1 GeV the 



box diagram contribution to the asymptotics is suppressed as it is seen from the behavior of the dashed line in Fig. [7] 
In the regime q 2 <§; lGeV 2 the MV asymptotics arises due to the pion-pole diagram discussed above (Fig. [3Jd). 
Indeed, we have 



Ait— Pole /*7(5 



Tr 


Q 2 X a 


Tr 







Gj 



1 - GpJpp(q 3 + qi) 



(64) 



,M4 



Nr. 



d 4 P 



(27T) 
d 4 k 

(2^r) 



4/ CP - ©) / (p + qi) Tr [T^Sip + Qi ) i l5 S (p - q 3 ) T^S (p)] 
J(k- qi )f(k + q 2 ) Tr [i To S (k - qi ) T pi S(k) T p2 S(k + q 2 )] 



The MV limit taken for the second integral corresponds to the asymptotics of the transition form factor for the pion 
of fixed virtuality considered earlier in [25|, l44j 



_c£k_ 

(2tt) 4 



/ (k - qi ) f (k + q 2 ) Tr [i l5 S (k - q x ) T m S (k) T^ 2 S(k + q 2 )] 



(65) 



2i 

<r 
2i 



^Suxsaq 6 N c 

q 2 



d 4 k 



f {k- qx) f(k + q 2 ) Tr [S (k - qi) i~f 5 S (k + q 2 ) 70-75] 



_- ^5„o- U (<l3) 
t v \Saq H3 rri r\ 21 

q 2 Tr [A 2 ] 



Thus the asymptotics of the LbL pion-pole diagram in the MV limit is 



— Pole 



P 5 = ^e PlP2 s p q*W^AT P3P + 



where 



AT p3P (q 3 ,q 4 )^q^U(q 2 )N c 



r Tr [T P3 S (p + q 4 ) i-y 5 S (p - q 3 ) T p S (p)] 



Gpf (p - q 3 ) f(p + q 4 ) 



(66) 



(67) 



(2tt) 4 l " " r "' J l-GpJ PP {q 3 ) 

This expression is similar to T^° l p c in ([59")) . If the renormalization of the axial vertex in (fBT)) coincided with AT^(k, k') 
in ([3"3")) then combining (|6"3"|) and ([6"6"| one would reproduce the MV asymptotics (fTB")) for all values of gf. However, 
the factor / w (53) has a strong dependence on the pion virtuality q\ and thus the required 1/q 2 dependence of the 
amplitude appears only from the pion-pole factor in ([67)) at q 2 <C 1 GeV 2 . 

The conclusions of this Section are: 1) We exactly reproduce the MV asymptotics (fT9|) in the regime q 2 3> 1 GeV 
from the quark box contribution to the LbL amplitude. 2) In the low-momentum region where effects of the infrared 
dynamics may be important we show how this asymptotics arises from the pion-pole triangle diagram near the position 
of the pole. 3) It might be possible that the MV asymptotics for arbitrary values of q 2 is reproduced within N^QM if 
we consider the contact terms contributing to the four-point amplitude, neglected in the present study. 4) It follows 
from the above results that we partially agree with MV when these authors state that the "pion-pole" contribution 
considered in [Tij is an attempt to describe the complete, on- and off-shell light-by-light scattering amplitude in the 
pseudoscalar isotriplet channel. 5) However, the NxQM model calculations do not support the factorization ansatz 
(ITU)) , with one form factor replaced by a constant, used by MV in order to simulate this complete amplitude 6 . 6) Finally, 

we agree with [4(| that in the MV model numerically one obtains = 7.97 • 10~ 10 instead of 7.65 • 10~ 10 quoted 



Moreover, the MV model does not fit the amplitude in other asymptotic perturbative QCD regimes: qf ~ q| ~ g| 3> ^qcd and 
1l ~ I3 I2 Aqcd' In NxQM when the photon virtualities are large the effective quark propagator connecting the hard photon 
vertices (7^) becomes the usual massless quark propagator. Thus, the asymptotics of the quark box diagram (Fig. 3a) is the same as 
in NxQM as well as in perturbative QCD for the massless quark considered in |l4|. In the ENJL model the quark box diagramm is 
suppressed at large photon virtualities by additional VMD form factors. 
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originally in fl4| . Presumably, the discrepancy has its origin in the numerical treatment of the large- momentum tails 
in the integrals. With the conformal variables, we treat these tails exactly. 

While the inclusion of the box diagram is crucial for the consistency of the NxQM, the numerical evaluation of 
its contribution to LbL is beyond the scope of this paper. We remark, that in the ENJL model [12] the result is 
2.1(0.3) • 10~ 10 , hence is a few times smaller than the dominant pion pole term. 



VIII. CONCLUSIONS 

In this paper we have analyzed the leading pion-pole contribution to a* - LbL in the nonlocal chiral quark model. 
The basic new element of our work is the inclusion of the full kinematic dependence of the pion-photon transition 
form factors, as it follows from the nonlocal chiral quark model. The dependence of the form factors on the pion 
virtuality decreases the result by about 15% compared to the case where this dependence is neglected. We have also 
demonstrated that the Melnikov-Vainshtein constraints, necessary for the consistency of the approach with QCD, are 
satisfied within the model when the quark box diagram is incorporated. Numerically, we quantitatively confirm the 
results obtained in other effective quark models. 

An important next step in the investigations of the light-by-light scattering within nonlocal quark model is to 
perform an extension to the so-called complete calculation (see EU HH ) , whch includes the scalar, axial- vector, 
and the 77, rf meson exchanges, as well as takes into account the quark and meson box diagrams. Due to contact terms 
arising in the nonlocal model, such a calculation is technically rather involved and will be presented elsewhere. 
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